clear all; clf;
map = [ ...
    0 6 4 6 2 6 4 6 0 ; ...
    6 5 6 5 6 5 6 5 6 ; ...
    4 6 3 6 4 6 3 6 4 ; ...
    6 5 6 5 6 5 6 5 6 ; ...
    2 6 4 6 1 6 4 6 2 ; ...
    6 5 6 5 6 5 6 5 6 ; ...
    4 6 3 6 4 6 3 6 4 ; ...
    6 5 6 5 6 5 6 5 6 ; ...
    0 6 4 6 2 6 4 6 0 ; ];
cm = colormap(jet(64));
cm2(1,:) = [.6 .6 .6];
cm2(2,:) = cm(55,:);
cm2(3,:) = cm(47,:);
cm2(4,:) = cm(40,:);
cm2(5,:) = cm(30,:);
cm2(6,:) = cm(20,:);
cm2(7,:) = [1 1 1];
colormap(cm2);
imagesc(map);
h = colorbar('ydir','reverse');
ytl = get(h,'yticklabel');
set(h,'ytick',.5:.85:6);
set(h,'yticklabel',ytl);
axis square